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1  Statement  of  Problem 

Animals  move  in  response  to  threats  from  predators,  a  change  in  the  environment,  or  to 
pursue  prey.  This  motion  is  generated  by  a  complex  interaction  of  the  neural  circuits  in  the 
spinal  cord  and  brain,  the  body,  the  muscles,  sensory  organs  and  the  environment.  With¬ 
out  sensory  input,  neural  circuits  are  capable  of  producing  movement  by  sending  signals  to 
motor  circuits  that  in  turn  contract  the  muscles.  Biomechanical  systems  of  body,  muscle 
and  environment  can  also  produce  motion  with  simple  activation  input.  However,  in  isola¬ 
tion,  neither  the  neural  circuit  or  the  biomechanical  system  can  produce  movement  efficient 
enough  to  support  chasing  prey  or  avoiding  predators.  Thus,  it  is  essential  that  the  neural 
circuit  receive  input  from  sensory  organs  that  are  affected  by  the  position  of  the  body,  the 
activation  of  the  muscle  and  the  environment.  Conversely,  efficient  movement  also  requires 
that  the  motor  system  receive  signals  from  the  neural  circuit.  Therefore,  a  central  ques¬ 
tion  in  modeling  integrated  neuromechanical  systems  is  how  to  close  the  loop  between  the 
neural  circuit  and  the  sensory  feedback  from  the  body,  muscle,  and  environment.  Under¬ 
standing  how  these  components  interact  at  a  systems  level  is  important  to  understanding 
the  complexities  in  vertebrate  locomotion. 

Vertebrate  locomotion  includes  walking,  swimming,  and  flying  and  each  of  these  be¬ 
haviors  involves  the  interaction  between  the  neural  circuit,  sensory  organs,  body,  muscle, 
and  environment.  Each  mode  of  locomotion  has  its  own  challenges  and  environments,  but 
the  feedback  loop  remains  similar.  We  have  closed  the  loop  between  the  neural  circuitry, 
body,  and  fluid  environment  for  swimming  locomotion,  where  the  lamprey  serves  as  a  model 
system 1'2'3,4,5’6’7’8. 

Our  philosophy  of  modeling  includes  a  hierarchy  of  models  with  varying  degrees  of  detail 
and  complexity.  Certainly  the  computational  power  to  simulate  an  all  neuron  simulation  with 
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full  muscle,  body  and  fluids  exists,  yet  the  results  may  not  yield  insight  into  the  mechanisms 
that  comprise  the  system.  Therefore,  we  have  chosen  to  start  with  a  tractable  system  in 
which  the  CPG  is  modeled  with  a  phase  model  and  the  plant  (body,  muscle  and  fluid)  is 
modified  from  McMillen  et  al.9.  It  is  crucial  to  understand  the  closed  loop  feedback  from 
the  simple  phase  model  perspective  before  incorporating  more  detailed  realistic  models  for 
CPG  or  the  fluid.  Looking  ahead  to  future  research  directions  that  would  build  on  the  work 
proposed  here,  the  next  step  in  the  modeling  hierarchy  is  to  combine  results  with  E.Tytell 
at  JHU,  who  is  using  Floquet  multiplier  theory  to  understand  which  eigenmodes  of  the  limit 
cycle  are  important  to  perturbations  in  the  environment.  Then,  the  combined  work  will  be 
integrated  into  a  two  or  three  dimensional  simulation  of  a  Navier-Stokes  fluid  and  clastic 
body3. 

2  Summary  of  results 

We  have  taken  two  approaches  to  understanding  the  role  of  sensory  feedback  in  locomotion: 
a  control-theory  approach  and  a  CPG-based  approach.  In  the  control  theory  approach  (§2.1), 
we  have  developed  a  general  method  of  computing  the  optimal  local  feedback  control  law 
given  a  model  of  the  plant  (the  mapping  from  muscle  activation  to  movement)  and  a  cost 
rate  function  that  penalizes  muscle  activation  and  rewards  speed  through  the  environment. 
The  first  step  in  this  method  is  to  the  compute  the  optimal  steady-state  swimming  pattern, 
which  we  have  done  for  a  simple  plant  model  of  lamprey  swimming  (§2.2).  In  our  CPG-based 
approach  (§2.3),  we  have  combined  a  simple  CPG  model  with  the  plant  model  of  §2.2  to 
produce  a  closed-loop  model  of  lamprey  swimming  in  which  a  the  CPG  produces  a  pattern 
of  muscle  activation  which  is  modified  by  sensory  feedback  related  to  body  curvature. 

2.1  Computation  of  optimal  feedback  for  the  control  of  locomotion 

We  have  a  developed  a  general  method  of  computing  the  optimal  feedback  control  law  for 
the  control  of  periodic  motion.  Consider  a  locomotion  control  problem  of  the  form 

x(t)  =  f(x(t),u(t)), 
cost  rate  g(x(t),u(t)), 

where  x(t)  is  the  state  vector,  u(t)  is  the  control  vector,  and  g(x,u )  punishes  the  use  of 
control  (e.g.,  muscle  activation)  and  rewards  speed  through  the  environment  in  the  desired 
direction. 

Let  u(t),  Xi(t),  x2(t)  be  the  open-loop  periodic  solution  that  minimizes  the  time  average 
of  g(u(t),x(t)).  Let  T  be  its  period  and  T  =  T/2tt.  Apply  the  change  of  coordinates  in  the 
vicinity  of  the  limit  cycle 

xi  =  X!  {T9), 
x2  =  x2  (T9)+r, 
u  =  u(T9)  +  u, 
g(x,u)  =  g{h{x,u)), 
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where  9  is  absolute  phase  (mod  2n),  r  describes  deviations  in  velocity,  u  describes  deviations 
in  the  control  signal,  x  =  (r,  9),  and  ( x,u )  =  h(x,u)  denotes  the  change  of  variables. 

In  the  new  variables  x  and  u,  the  system  has  the  form 

f  \  =  it*  =  (  Ava (0)r  +  A01(6)u  +  A20(9)r2  +  Au(9)ru  +  A02{9)u2  \  n(  , 

9  )  nx*U)  \  u  +  B10(9)r  +  B01(9)u  +  B20(9)r 2  +  Bu(9)ru  +  B02(9)u 2  )  + 


g(x,u)  =  g00(9)  +  gw{9)r  +  g0i{9)u 

+  g2o{0)r2  +  gu(9)ru  +  g02(9)u2  +  0(3), 

where  u  =  1/T  and  we  have  the  subtracted  the  mean  from  goo(O)- 

For  any  initial  state  (r(0)),  ^(0)),  we  want  to  choose  u(t)  for  0  <  t  <  tf,  where  tf  is  a  time 
far  in  the  future,  to  minimize 

l/(r(O),0(O))  =  D(9(tf))  +  [  g(r(t),9(t),u(t))dt. 


where 


D(9)  = - /  g00{(p)d(j) 

u  Jo 


so  that  V  converges  in  the  limit  as  tf  — >  oo. 

The  solution  is  given  by  the  Hamilton- Jacobi-Bellman  (HJB)  equation 

0  =  min[g(r,  9,  u )  +  W (r,  9)  ■  /(r,  6*,  u)] 

U 

Using  a  perturbation  analysis,  we  can  approximate  the  solution  of  the  HJB  equation: 

V(r,9 )  =  Vq(0)  +  V\(9)r  +  rTV2(9)r  +  0(|  |r|  |3), 
u  =  Ci(9)r  +  0(r2), 


so  that  the  solution  of  optimal  feedback  control  problem  in  the  vicinity  of  the  limit  cycle  is 

u  =  C0{9)  +  CMr  +  0{r2) 

where  Cq(9)  =  u(T9). 

The  quantity  F(r,9,u)  =  g(r,9,u )  +  W(r,9 )  •  f(r,9,u )  to  be  minimized  in  the  HJB 
equation  can  be  written  up  to  second  order  as 

F(r,  9,  u)  =  F00(9)  +  F10{9)r  +  F01(9)u  +  rTF20(9)r  +  rTFu(9)u  +  uTF02(9)u 
=  goo(0)  +  gw(9)r  +  goi(9)u  +  rTg20(9)r  +  rTgn(0)u  +  uTg02(9)u 
+  (hi (9)  +  2rTV2(9)) (A10(9)r  +  A01(9)u  +  A20(9)[r ,  r]  +  Au(9)[r ,  u]  +  A02(9)[u,  u\) 

+  (V'(9)  +  rTVl(9)T  +  rTV'(9)r) 

(u>  +  B10(9)r  +  BOi(0)u  +  rTB20(9)r  +  rTBn(0)u  +  uTB02(9)u), 
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where 


Foo(0)  —  ^oo(^)+c oVq(9), 

Fw(6)  =  gw(0)  +  Vi(9)Aw(9)  +  Vq(9)Biq(9)  +  c oV[{9), 

Fm{9)  =  goi{9)  +  VmAoi{9)  +  V^9)Bn{9), 

F2q{9)  =  cj2o(9)  +  (hi  o  A2O)(0)  +  2V2 (9)A10(9)  +  Vq(9)B20(9)  +  V^(0)TB1O(9)  +  wV2(9), 
Fn(0)  =  gn(0)  +  (V.o  An)(0)  +2V2(0)AO1(0)  +  V'(0)BU(0)  +  VI(0)TBO1(0),  (1) 

F02(9)  =  go2(9)  +  (V1oAO2)(9)  +  V'(9)Bo2(0).  (2) 

To  find  the  u  that  minimizes  F(r,0,u )  for  each  r  and  u ,  we  solve 

0  =  VuF(r,9,u)  =  Fol(0f  +  Fll{0)Tr  +  2  FO2(0)u 

for  u.  For  r  =  0,  the  solution  must  be  u  =  0,  since  this  corresponds  to  the  solution  of  the 
open-loop  optimal-control  problem.  Therefore,  F0l(9)  =  0.  Then  the  solution  to  (2.1)  is 

u  =  ^FvW-'FuWfr.  (3) 

Therefore, 

minF(r,  0,  u)  =  Foo(0)  +  Flo(0)r  +  rTF20(9)r  -  lrTF11{0)FO2(0y1Fll(0)Tr 

U 

+  yF11(0)FO2(0)-1Fu(0)Tr 

=  Foo(0 )  +  F10(9)r  +  rT  (F20(9)  -  lF11(0)FO2(9y1F11(9)T)  r. 

To  solve  the  HJB  equation,  we  want  this  minimum  to  be  0  for  all  r  and  9.  Therefore, 

F00(9)  =  0, 

Fi  o(0)  =  0, 

F20(9)  -  1^11(0)^02 (9)-1Fu(0)T  =  0. 

We  write  F2 o(0)  and  Fn(9)  in  the  form 

F2O(0)  =  G2O(0)  +  2V2(0)Alo(0)  +  wV'(0), 

Fn(0)  =  G11(0)  +  2V2(0)AO1(0), 

where 

G2o(0 )  =  g2o (9)  +  (h'i  °  A2o)(9)  +  Vo'(0)B2O(0)  +  V1'(9)B10(9), 

Gn(0)  =  gn(9)  +  (V1oAu)(9)  +  V'(0)Bu(0)  +  V;(9)Bol(9). 

Then  V2(9)  is  the  solution  to 

G2o(d)  +  2V2(0)AiO(0)  +  wV2{0) 

-  i[Gn(0)  +  2V2{9)Am{9)]F02{9)-l[Gll{9)  +  2V2(0)Ao1(0)}t  =  0. 

Once  we  have  solved  for  V0(9),  V\{9)  and  V2(0),  we  can  use  (1)  and  (2)  to  solve  for  Fn(9) 
and  F02(9)  and  then  use  (3)  to  obtain  the  optimal  feedback  control  law  in  the  vicinity  of  the 
limit  cycle. 
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2.2  Computation  of  optimal  steady-state  swimming  for  a  simple 
model  of  lamprey  swimming 

To  apply  the  method  of  §2.1  to  a  plant  model  of  lamprey  swimming,  the  first  step  is  to 
find  the  optimal  steady-state  swimming  solution.  We  have  done  this  for  a  simple  model 
of  lamprey  swimming,  as  described  in  this  subsection.  (We  are  currently  in  the  process  of 
implementing  the  feedback  step  of  the  method  of  §2.1  for  this  model.) 

The  plant  model  we  used  is  based  on  Ekeberg  (1993) 1  and  McMillen  &  Holmes  (2006) 9 
with  further  simplifications.  The  body  is  modeled  as  a  chain  of  n  segments  with  n  —  1  links. 
For  simplicity,  we  assume  that  all  segments  are  identical  (in  contrast  to  the  lamprey  body 
which  tapers  from  head  to  tail).  The  system  has  n  +  2  mechanical  degrees  of  freedom:  the 
2-d  position  of  the  head  and  the  angles  of  each  of  the  n  segments.  Therefore,  the  state 
is  described  by  a  vector  x(t )  of  dimension  2 n  +  4.  We  assume  that  when  a  segment  moves 
through  the  water  there  are  resistive  water  forces  in  the  directions  perpendicular  and  parallel 
to  the  segment’s  orientation  that  depend  on  the  segment’s  velocity  components  in  these  two 
directions. 

We  assume  that  the  control  vector  u(t)  consists  of  the  torques  produced  by  muscles  at 
each  of  the  n  —  1  links.  Therefore,  we  have  a  plant  model  of  the  form 

x(t)  =  f(x(t),u(t)), 

where  x{t)  has  dimension  2 n  +  4  and  u{t)  has  dimension  n  —  1. 

We  assume  a  simple  cost  rate  that  penalizes  mean  squared  muscles  torques  and  rewards 
velocity  of  the  center  of  mass  ucom(^)  hi  the  leftward  direction: 

^  n—  1 

cost  rate  g(x(t),  u(t))  =  -  cucom(^), 

H  1  . 

J=1 

where  the  parameter  c  specifies  the  relative  importance  of  swimming  fast  versus  minimizing 
muscle  activation. 

To  optimize  steady-state  swimming  we  minimize  the  average  cost  rate  over  one  cycle  of 
swimming  by  adjusting  the  following  set  of  Qn  —  2  parameters:  the  swimming  frequency 
fs  =  1/T;  the  initial  2-d  velocity  of  the  head  and  time  derivative  of  the  first  segment  angle; 
and  the  Fourier  coefficients  for  modes  1,  3  and  5  describing  link  angles  as  a  function  of  time. 
We  assume  left-right  spatio-temporal  symmetry,  so  we  only  need  the  odd-numbered  modes. 
We  assume  that  that  at  t  —  0,  the  2-d  head  position  and  first  segment  angle  are  zero.  We 
have  the  equality  constraints  that  the  2-d  velocity  of  the  head  and  time  derivative  of  the 
first  segment  angle  must  be  the  same  at  the  beginning  and  end  of  the  cycle.  To  reduce  the 
number  of  parameters,  here  we  set  the  number  of  links  n  to  be  6.  We  solve  the  constrained 
optimization  problem  using  the  Lagrange  method  as  implemented  by  the  fmincon  function 
in  Matlab.  The  average  cost  rate  is  computed  by  numerically  integrating  the  plant  system 
of  ODEs  over  one  cycle. 

Figure  1  shows  results  of  optimizing  swimming  for  different  rewards  c  for  swimming 
speed.  As  the  reward  c  increases,  the  lamprey  swims  faster  (Fig.  1A).  As  in  the  lamprey,  as 
swimming  speed  in  the  model  increases,  swimming  frequency  increases  (Fig.  IB)  and  phase 
lags  among  muscle  torques  and  body  curvature  at  different  locations  along  the  body  remain 
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roughly  constant  (Fig.  1C).  Also  note  that  as  in  the  lamprey,  the  wave  of  muscle  activity 
moves  more  faster  along  the  body  than  the  wave  of  body  curvature.  As  expected  muscle 
torques  also  increase  with  increasing  swimming  speed  (Fig.  ID). 

2.3  A  closed-loop  CPG-based  model  of  lamprey  swimming 

In  this  subsection  we  describe  our  CPG-based  approach  to  understanding  closed-loop  lo- 
comtion.  Vertebrate  locomotion  central  pattern  generators  (CPG)  are  known  to  produce  an 
electrical  signal  that  is  projected  to  the  motornenrons,  which  in  turn  contract  the  muscle 
and  propel  the  animal  through  the  environment.  Sensory  input  to  the  vertebrate  locomo¬ 
tion  CPG  is  known  to  modulate  the  signal  of  the  CPG.  In  the  model  system  that  we  are 
studying,  the  lamprey  have  several  forms  of  sensory  inputs.  The  lamprey  lateral  line  system 
transmits  information  regarding  the  fluid  flow  around  the  body  as  a  direct  input  through 
the  lateral  line  nerve  to  the  brain.  Additionally,  there  are  edge  cells  that  sit  on  the  margin 
of  the  spinal  cord  that  transmit  information  directly  to  the  locomotion  CPG.  Edge  cells  are 
stretch  receptor,  and  essentially  measure  the  curvature  of  the  body.  In  our  model  of  closed 
loop  swimming,  we  focus  on  this  type  of  sensory  input  and  show  that  sensory  input  to  the 
lamprey  locomotion  CPG  modulate  the  motion  of  the  animal  though  a  closed  feedback  loop. 

The  vertebrate,  and  lamprey  in  particular,  central  pattern  generator  of  locomotion  is 
commonly  represented  by  a  chain  of  coupled  oscillators  with  neural  connections  between  the 
oscillators.  Each  individual  oscillator  can  be  represented  by  models  with  varying  biological 
detail6,10’1,8’11.  We  consider  the  simplest  model  of  a  chain  of  n  coupled  phase  oscillators: 

n 

Oi  =  co  +  ^2  ai~i  sin(0j  _  Oi  ~  A-j)  +  A i  sin(6>“  -  0*  -  -0D,  (4) 

3= 1 

where  9t  for  i  =  1 , ,  n  are  the  absolute  phases  (mod  2n)  of  oscillators  in  the  chain  (?'  =  1 
for  the  oscillator  nearest  the  head)  and  the  parameter  u  is  the  uncoupled  angular  frequency 
of  the  oscillators  in  the  chain,  which  we  assume  to  be  the  same  for  all  oscillators,  u  =  1  Hz. 

Each  oscillator  in  the  central  pattern  generator  receives  two  types  of  input.  The  first  type 
of  input,  which  we  will  call  internal  inputs ,  are  from  other  oscillators  in  the  central  pattern 
generator.  This  describes  the  coupling  between  the  oscillator  in  the  locomotion  CPG.  The 
model  assumes  that  coupling  is  translation  invariant ,  meaning  that  the  properties  of  the 
connection  from  oscillator  j  to  oscillator  i  only  depend  on  the  relative  position  k  =  i  —  j: 
ak  is  the  strength  of  the  connection  and  ipk  is  its  preferred  phase.  Here,  k  >  0  for  de¬ 
scending  (toward  the  tail,  or  caudal)  connections  and  k  <  0  for  ascending  (toward  the  head, 
or  rostral)  connections.  (Note  the  conflict  between  the  anatomical  terms  “descending”  and 
“ascending”  and  our  indexing  scheme  in  which  i  increases  from  head  to  tail.)  There  is  no 
connection  between  an  oscillator  and  itself,  so  a0  =  0.  Supported  by  numerical  evidence12,13, 
we  assume  all-to-all  coupling  with  coupling  strength  decaying  exponentially  with  distance 
|r|:  for  r  >  0  (descending  connections),  ar  =  Ade~^r^Xd]  for  r  <  0  (ascending  connections), 
ar  =  Aae_lrl'/A“;  for  r  =  0  (intrasegmental  connections),  ar  =  1,  where  Ad,  Xd  and  Aa,  Xa 
are  the  amplitudes  and  length  constants  for  descending  and  ascending  connections,  respec¬ 
tively.  Experimental  evidence  suggests  non-uniform  coupling  asymmetry,  specifically,  that 
short  ascending  connections  are  stronger  than  short  descending  connections,  whereas  long 


7 


descending  connections  are  stronger  than  long  ascending  connections,  which  are  reflected  in 
the  parameters  Ad,  A  a  and  Aa,  Xa. 

The  second  type  of  input  into  the  central  pattern  generator  is  external  inputs  from  the 
edge  cells.  This  input  provides  the  sensory  feedback  information  into  the  central  pattern 
generator.  Experimentally,  it  has  been  shown  that  edge  cells  project  directly  onto  the  central 
pattern  generator,  and  that  it  responds  to  both  stretch  and  rate  of  stretch  of  the  body,  but 
the  exact  form  of  the  feedback  function  has  not  been  experimentally  determined.  We  model 
the  input  from  the  edge  cells  using  curvature  information  from  the  body  equations.  We 
tried  two  forms  of  feedback  functions:  one  that  was  proportional  to  the  body  curvature,  and 
one  that  was  proportional  to  the  square  of  the  body  curvature.  We  determined  values  of 
the  constant  of  proportionality,  which  we  will  call  the  strength  of  feedback  and  denote  by 
af.  for  which  the  animal  reaches  steady  state  swimming.  The  feedback  function  that  was 
proportional  to  the  body  curvature  produced  unrealistic  ranges  of  feedback  strength  that 
produced  steady  state  swimming.  The  ranges  in  feedback  strength  to  achieve  steady  state 
swimming  were  on  the  order  of  10~6,  which  essentially  represents  no  feedback  from  the  edge 
cells  onto  the  central  pattern  generator.  The  feedback  function  that  was  proportional  to 
the  square  of  body  curvature  produced  more  realistic  results  in  terms  of  values  of  feedback 
strength  that  produced  steady  state  swimming,  and  these  results  are  shown  below. 

For  our  simulations,  we  used  a  lamprey  with  (n  =  10)  segments,  a  Young’s  modulus  of 
0.7,  and  all-to-all  coupling  between  the  CPG  oscillators  with  a  slight  non-uniform  coupling 
asymmetry  with  Ad  =  10.0  and  Aa  =  10.1.  The  simulation  was  performed  in  Matlab  using 
the  ode23s  package.  Figure  2a  illustates  the  phase  lags  relative  to  the  first  oscillator  for 
oscillators  2,4,6,  and  8  of  the  body,  in  red,  blue,  green,  cyan,  respectively.  The  effect  of 
the  feedback  is  to  slightly  increase  the  phase  lags  between  the  oscillators,  with  oscillator  8 
showing  the  largest  increase,  as  expected  since  in  incorportates  the  previous  six  increases. 
Figure  2b  shows  that  the  period  is  constant,  as  a  function  of  body  position,  as  expected  in 
steady  state  swimming,  but  the  period  is  different  for  each  value  of  the  feedback  strength. 
In  particular,  the  period  decreases  as  the  feedback  strength  increases.  Figure  2c  shows 
the  maximum  curvature  as  a  function  of  body  position  for  varying  values  of  the  feedback 
strength.  As  the  feedback  strength  increases,  the  maximum  curvature  decreases  for  each 
body  position,  although  the  difference  between  the  maximum  curvatures  at  points  along  the 
body  for  a  given  feedback  strength  seems  to  increase  as  feedback  strength  increases.  The 
final  graph,  figure  2d  shows  swimming  speed  as  a  function  of  feedback  strength.  There  is 
a  sharp  rise  in  swimming  speed  for  small  feedback  strength.  The  sharp  increase  is  followed 
by  a  gradual  decrease  in  swimming  speed  as  feedback  strength  increases.  We  note  here  that 
swimming  speeds  on  the  order  of  10-40  cm/sec  are  not  realistic  for  a  swimming  lamprey, 
which  may  indicate  that  the  feedback  function  is  more  complicated  than  simply  squared 
curvature. 


Phase  lag  relative  to  the  first  oscillator 


Feedback  Strength  Position  Along  Body 


(a)  Phaselags  Relative  to  the  first  Oscillator 


(b)  Period 


(c)  Curvature 


(d)  Swimming  Speed 


Figure  2:  Closed-loop  swimming  for  various  levels  of  feedback.  This  figure  illustrates  the 
effect  of  feedback  on  phase  lags  relative  to  the  first  oscillator  in  figure  2a,  period  in  figure  2b, 
the  curvature  as  a  function  of  body  position  in  figure  2c  and  velocity  of  swimming  in  figure  2d. 
In  figure  2a,  the  phaselags  relative  to  the  first  oscillator  are  shown  as  a  function  of  feedback 
strength  otf,  in  a  range  of  values  that  produce  steady  state  swimming.  Relative  phaselags 
are  shown  for  oscillator  2  in  red,  oscillator  4  in  blue,  oscillator  6  in  green  and  oscillator  8  in 
cyan.  Figure  2b  shows  that  the  period  is  a  constant  function  of  body  position,  as  expected  in 
steady  state  swimming,  but  that  period  decreases  as  feedback  strength  increases.  Figure  2c 
plots  maximum  curvature  as  a  function  of  body  positions  for  varying  values  of  feedback 
strength.  The  body  curvature  decreases  as  a  function  of  feedback  strength.  The  final  graph 
figure  2d  shows  the  effect  of  feedback  strength  on  swimming  speed. 
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